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The steady two-dimensional laminar stream of a gooey incompressible electrically directing liquid over an 
exponentially extending sheet in the presence of magnetic field with viscous dissipation, heat generation, chemical reaction 
and radiative heat flux is studied. By suitable similarity transformation, the overseeing limit layer comparisons are changed 
to normal differential mathematical statements and illuminated numerically by utilizing fourth request Runge-Kutta 
technique with shooting method. The impacts of different overseeing parameters on the speed, temperature, focus, 
skin-erosion coefficient, Nusselt number and Sherwood number are registered and examined in subtle element. 
Comparisons with previously published work are performed and are found to be in a good agreement. 
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In the most recent couple of decades, liquid stream with warmth and mass exchange on a constantly extending 
sheet has pulled in extensive consideration due to its numerous applications in businesses, designing and assembling 
methods. Illustrations of these applications incorporate the glass-fiber creation, wire drawing, paper generation, plastic 
sheets, metal and polymer handling businesses, hot moving and constant throwing of metals and turning of filaments. 
The kinematics of extending and the synchronous warming or cooling amid such methods assume an imperative part on the 
structure and nature of the last item. Numerous scientists enlivened by Sakiadis [1,2] who launched the limit layer conduct 
examined the extending stream issue in different viewpoints. Augmentation to that, a careful arrangement was given by 
Crane [3] for a limit layer stream brought on by extending surface. 

Gupta [4], Carragher and Crane [5], Dutta et al. [6] mulled over the warmth move in the stream over an extending 
surface considering distinctive parts of the issue. Magyari and Keller [7] watched that the investigation of limit layers on 
an exponentially extending consistent surface with an exponential temperature circulation. Sanjayanand and Khan [8] 
studied the heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet. They found 
that the viscoelastic parameter enhances the thermal boundary layer thickness. The effect of viscous dissipation on the 
mixed convection heat transfer from an exponentially stretching surface was studied by Partha et al. [9]. They observed a 
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rapid growth in the non-dimensional skin friction coefficient with the mixed convection parameter. 

The attributes coveted of the last item in an expulsion methodology rely on upon the rate of extending and 
cooling. Subsequently, it is critical to have a controlled cooling environment where the stream over the extending sheet can 
be directed by outside offices like an attractive field. An exponential variety of an attractive field is utilized, among 
different applications, to focus the diamagnetic weakness of plasma. Pavlov [10] considered the magnetohydrodynamic 
stream of an incompressible gooey liquid over a straightly extending surface. Sarpakaya [11] extended Pavlov's work to 
non-Newtonian liquids. Ensuing studies by Andersson [12], Lawrence and Rao [13], Abel et al.[14], Cortell [15] concerned 
the magnetohydrodynamic stream of viscoelastic fluids over an extending sheet. 

The majority of the prior work dismissed radiation impacts. In the event that the polymer expulsion methodology 
is put in a thermally controlled environment, radiation could get to be vital. Numerous analysts have considered the impact 
of warm radiation on streams over extending sheets. The impact of warm radiation on the limit layer stream because of an 
exponentially extending sheet is examined by Sajid and Hayat [16]. Studies by Raptis [17], Raptis and Perdikis [18] 
address the impact of radiation in different circumstances. The impacts of radiation on hydromagnetic limit layer stream of 
a consistently extending surface have pulled in extensive consideration as of late because of its various applications in 
industry. Kameswaran et al.[19] watched that radiation impacts on MHD Newtonian fluid stream because of an 
exponential extending sheet. Seini and Makinde [20] mulled over impacts of radiation and substance response on MHD 
limit layer stream over exponential extending surface. Siddheshwar and Mahabaleswar [21] concentrated on the impacts of 
radiation and warmth source on MHD stream of a viscoelastic fluid and warmth exchange over an extending sheet. Bidin 
and Nazar [22] contemplated the impacts of numerical arrangement of the limit layer stream over an exponentially 
extending sheet with warm radiation. Elbashbeshy and Dimian [23] investigated limit layer stream in the vicinity of 
radiation impact and warmth exchange over the wedge with a gooey coefficient. Warm radiation consequences for 
hydro-attractive stream because of an exponentially extending sheet were mulled over by Reddy and Reddy [24]. Raptis et 
al. [25] examined the impact of warm radiation on the magnetohydrodynamic stream of a gooey liquid past 
semi-unbounded stationary plate and Hayat et al. [26] amplified the examination for the second grade liquid. Jat and Gopi 
Chand [27] found that the effects of dissipation and radiation on MHD flow and heat transfer over an exponentially 
stretching sheet. Radiation effects on MHD flow past an exponentially accelerated isothermal vertical plate with uniform 
mass diffusion in the presence of a heat source was studied by Reddy et al. [28]. They observed that the speed abatements 
with an increment in the attractive parameter because of a resistive drag power which has a tendency to oppose the liquid 
stream and along these lines decreases the speed. The limit layer thickness was additionally found to decline with an 
increment in the attractive parameter. 

Notwithstanding an attractive field and warm radiation, one needs to consider the gooey dissemination impacts 
because of frictional warming between liquid layers. The impact of thick dissemination in characteristic convection 
techniques has been examined by Gebhart [29] and Gebhart and Mollendorf [30]. They watched that the impact of thick 
dissemination is dominating in overwhelming regular convection and blended convection forms. They additionally 
demonstrated the presence of a similitude answer for the outer stream over an endless vertical surface with an exponential 
variety of surface temperature. 

The investigation of warmth era or assimilation in moving liquids is critical in issues managing compound 
responses and those concerned with separating liquids. Heat era impacts may adjust the temperature dissemination and this 
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thusly can influence the molecule testimony rate in atomic reactors, electronic chips and semi conductor wafers. Albeit 
definite displaying of inward warmth era or ingestion on is truly troublesome, some basic scientific models can be utilized 
to express its general conduct for most physical circumstances. Vajravelu and Hadjinicalaou [31] studied the heat transfer 
characteristics over a stretching surface with viscous dissipation in the presence of internal heat generation or absorption. 
Mohammed Ibrahim and Bhaskar Reddy [32] noticed that effects of radiation and mass transfer effects on MHD flow 
along a stretching surface in presence of viscous dissipation and heat generation. Pavithra and Gireesha [33] noticed that 
the effect of heat generation on dusty fluid flow over an exponentially stretching sheet with viscous dissipation. 
Mohammed Ibrahim [34] found that radiation effects on mass transfer flow passing through a highly porous medium in 
presence of heat generation and chemical reaction. Mohammed Ibrahim [35] investigated effects of chemical reaction and 
radiation on MHD free convection flow along a stretching surface in presence of dissipation and heat generation. 

In the present article, we investigate the effects of various physical and fluid parameters such as the magnetic 
parameter, radiation parameter, viscous dissipation parameter, heat generation parameter and chemical reaction parameter 
on the flow, heat and mass transfer characteristics of an exponentially stretching sheet. The force, vitality and fixation 
mathematical statements are coupled and nonlinear. By utilizing suitable comparability variables, these comparisons are 
coupled and nonlinear. By utilizing suitable comparability variable, these comparisons are changed over into coupled 
customary differential mathematical statements and are illuminated numerically by utilizing shooting procedure with the 
forward request Range-Kutta strategy. 

Formulation of the Problem 




Figure 1: Schematics of the Problem 

Consider an enduring two dimensional laminar stream of a thick incompressible electrically leading liquid over a 
continuous exponentially stretching surface. The x - axis is taken along the stretching surface in the direction of motion 
and y- axis is perpendicular to it. The sheet velocity is assumed to vary as an exponential function of the distance x from 

T C 

the slit. The temperature and concentration far away from the fluid are assumed to be °° and °° respectively as shown in 
Figure 1. The sheet-ambient temperature and concentration differences are also assumed to be exponential functions of the 
distance x from the slit. A variable magnetic field of strength B(x) is applied normally to the sheet. Under the usual 
boundary layer approximation, subject to radiation, viscous dissipation, heat generation and chemical reaction effects, the 
equations governing the momentum, heat and mass transports can be written as 

Continuity equation 
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du dv 
— + — = 0 
dx dy 



Momentum equation 



du dv d u <jB a 

u—+v— = v — —u 

ox dy dy p 



(1) 



(2) 



du 



Energy equation 

dT dT k d 2 T 1 dq r V 
u — + v — = — + — 

dx dy pc p dy pc p dy c p [ dy 

Species equation 
dC dC „d 2 C *, v 

" & +v * = D W~ Kr (c_c - ) 



2 



+ ?B_ u 2 + Q^( T _ T J 



(3) 



(4) 



where w and v are the velocity components in the x, y directions respectively, ^ is the kinematic viscosity, P is 
the density, & is the electrical conductivity of the fluid, T is the temperature, C is the concentration, k is the thermal 

conductivity, c p is the specific heat at constant pressure, g r is the radiative heat flux, ^° is the heat generation coefficient, 
D is the species diffusivity, Kr j s the reaction rate parameter. 

The boundary conditions for the velocity, temperature and concentration fields are 

x 2x 2x 

u = U w =U 0 e L v = 0 T = T w =T„+T 0 e T C = C w = C„ + C 0 e^ 



u^O T-*T„ C-+C„ v-> 



at y = 0 

as (5) 

oo T C 

Here the subscripts ' refer to the surface and ambient conditions respectively, 0 ' 0 are positive constants, 

^° is the characteristic velocity, L is the characteristic length and Q is the constant. 

To facilitate a similarity solution, the magnetic field B(x) is assumed to be of the form 

B(x) = B 0 e~ L (6) 

where B 0 is a constant. It is also assumed that the fluid is weakly electrically conducting so that the induced 

Q 

magnetic field is negligible. Following Rosseland's approximation, the radiative heat flux r is modeled as 
\o" dT A 

3k* dy (7) 
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where & is the Stefan-Boltzman constant, is the mean absorption coefficient. 

This estimate is substantial at focuses optically a long way from the limit surface and it is useful for serious 
ingestion, which is for an optically thick limit layer. It is expected that the temperature contrast inside the stream such that 

>^i... may be communicated as a straight capacity of temperature. Subsequently, growing 1 by Taylor series 
T 

about °° and disregarding higher-order terms 
T 4 = \fj - 3T* 

dg r _ 16o"X 3 d 2 T 

We have 3 ? 3 ^ V (8) 

Continuity equation (1) is satisfied by introducing a stream function ^ such that 



u -—— v = — 

oy dx 



The following similarity variables are used: 



(9) 



u = U 0 e L f(rj) V V 2L e 2L {fW + Vf '(*?)) 



T = T ao +T 0 e L 0(7]) C = C„ + C 0 e L m * = V^T ye* 



u 0 



(10) 



Where ^ is the similarity variable, f^l^ is the dimensionless stream function, is the dimensionless 

temperature, and [ s the dimensionless concentration. 

On using equations (6), (8) and (10), equations (2) - (5) are transformed to: 

r-2f' 2 + ff"-Mf = 0 
f 4 > 



1 + -R 

3 



9" + Pr [/(?'- f'S + Ec(f) 2 +MEc(f) 2 + Q6~\ = 0 

(12) 



f + Scf(f>' - Scf(/) - ScKr(/> = 0 



(13) 



/(0) = 0, /'(0)=1 ; /'(°°)->0, (14) 
0(0) = 1, 0(°o)^O, (15) 

m=i 0(°°)^>o. (16) 

The non-dimensional constants showing up in mathematical statements (11) - (13) are the magnetic parameter M, 
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the radiation parameter R, the Prandtl number Pr, the Eckert number Ec, Q is the heat generation parameter, Sc is the 
Schmidth number and Kr is the chemical reaction parameter individually characterized 



M = 2aB ° L 

pu 0 



R 



4<rX 3 Pr = ^JL Ec = ^ 



k*k 



2Q 0 L 2LKr * 

— Sc = — Kr = — 

D U 0 



U 0 e L 



Skin Friction, Heat and Mass Transfer Coefficients 

The parameters of designing enthusiasm for warmth and mass transport issues are the skin erosion coefficient, the 
neighborhood Nusselt number Nu, and the nearby Sherwood number Sh. These parameters separately describe the surface 
drag, divider warmth and mass exchange rates. 

The shearing stress at the surface of the wall w is given by 
ou 



juU 0 Re 



3x 



L \ 2 e 2L f"(0) 



(17) 



where ^ is the coefficient of viscosity and ^ is the Reynolds number. The skin friction coefficient is 



defined as 



2^ 



(18) 



and using equation (17) in equation (18), we obtain 



= -/*(0) 



(19) 



The heat transfer rate at the surface flux at the wall is given by 



4w 



_-k(T w -T„) Re 



2 e 2L 6\Q) 



(20) 



where k is the thermal conductivity of the fluid. The Nusselt number is defined as 
Nu = *^- 

Using Equation (20) in Equation (21), the dimensionless wall heat transfer rate is obtained as follows: 



(21) 
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Nu =-e\o) 

/Re/ 

2 ■ (22) 

The mass flux at the surface of the wall is given by 

_-D(C w -CJ /Re" , 



r dC^ 



2 e 2L 0'(O) 



(23) 



The Sherwood number is defined as 

DC W -C„ (24) 

Using (23) in (24), the dimensionless wall mass transfer rate is obtained as 

Sh ■ = -*'«» 



Re, 



2 (25) 
In Equations (19), (22) and (25), Re represents the local Reynolds number and it is defined as 

V 

Numerical Procedure 

The set of nonlinear ordinary differential equations (11), (12), and (13) with boundary conditions (14) - (16) were 
solved numerically using Runge - Kutta fourth order algorithm with a systematic guessing of f , ^ and ^ 
by the shooting technique until the boundary conditions at infinity are satisfied. The step size = 0.001 is used while 

lxl(T 5 

obtaining the numerical solution and accuracy up to the fifth decimal place i.e. ' which is very sufficient for 

71 — > °° V 

convergence. In this method, we choose suitable finite values of ' , say '°° , which depend on the values of the 

parameter used. The computations were done by a program which uses a symbolic and computational computer language 
in Mathematica. 

RESULTS AND DISCUSSIONS 

To dissect the outcomes, numerical calculation has been done utilizing the system depicted as a part of the past 
passage for different in representing parameters, in particular, magnetic field parameter M, Prandlt number Pr, radiation 
parameter R, heat generation parameter Q, Eckert number Ec, Schmidt number Sc, chemical reaction parameter Kr. In the 
present study taking after default parameter qualities are received for processings:M = 1.0, Pr = 0.71, R = 0.5, Q = 0.1, Ec 
= 0.1, Sc = 0.6, Kr = 0.5. All graphs therefore correspond to these values unless specifically indicated on the appropriate 
graph. 

Figure 2 shows the variation of the velocity profile against the magnetic parameter. We notice that the effect of 
the magnetic parameter is to reduce the velocity of the fluid in the boundary layer region. This is due to an increase in the 
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Lorentz force, similar to Darcy's drag observed in the case of flow through a porous medium. This adverse force is 
responsible for slowing down the motion of the fluid in the boundary layer region. 

The variation of the temperature distribution with the magnetic parameter is shown in Figure 3. The thermal 
boundary layer thickness increases with increasing values of the magnetic parameter. The opposing force introduced in the 
form of the Lorentz drag contributes in increasing the frictional heating between the fluid layers, and hence energy is 
released in the form of heat. This results in thickening of the thermal boundary layer. 

The effect of the magnetic parameter on the concentration profile is shown in Figure 4. It is observed that 
increases in the values in M result in thickening of the species boundary layer. 

Figure 5 illustrates the temperature profile for different values of the Prandtl number Pr. The Prandtl number 
defines the ratio of momentum diffusivity to thermal diffusivity. The numerical results show the effect of increasing values 
of Prandtl number results in decreasing temperature. The reason is that smaller values of Prandtl number Pr are equivalent 
to increasing the thermal conductivities, and therefore heat is able to diffuse away from the heated plate more rapidly than 
for higher values of Prandtl number Pr. Hence in the case of smaller Prandtl number Pr as the boundary layer is thicker 
and the rate of heat transfer is reduced. 

The influence of the thermal radiation parameter R on temperature is shown in Figure 6. It is clear that thermal 
radiation enhances the temperature in the boundary layer region. Thus radiation should be kept at its minimum in order to 
facilitate better cooling environment. The radiation parameter R defines the relative contribution of conduction heat 
transfer to thermal radiation transfer. 

The effect of the Eckert number Ec on heat transfer is shown in Figure 7. It is clear that the temperature in the 
boundary layer region increases with an increase in the viscous dissipation parameter. 

Figure 8 shows the influence of the heat generation parameter Q on the temperature profile within the thermal 
boundary layer. From the Figure 8 it is observed that the temperature increases with an increase in the heat generation 
parameter. 

Figures 9-10 depict chemical species concentration profiles against co-ordinate ^ for varying values physical 
parameters in the boundary layer. The species concentration is highest at the plate surface and decreases to zero far away 
from the plate satisfying the boundary condition. From these figures, it is noteworthy that the concentration boundary layer 
thickness decreases with an increase in Schmidt number and chemical reaction parameter. 

We likewise take note of that since the vitality mathematical statement is incompletely decoupled from the force 
and species preservation comparisons, the parameters influencing the vitality comparison, specifically, the Prandtl number, 
the thermal radiation parameter, heat generation parameter and the Eckert number, don't change velocity and concentration 
profiles. Table 1 shows the comparison of Kameswaran et al. [19] work with the present work for Pr = R = Ec = Q = Sc = 
Kr = 0 and it note worthy that there is a good agreement. 

Table 2 indicates the values of skin-friction coefficient, the wall temperature gradient and the wall concentration 

gradient in terms of f ® ^ and ^ reS p ec ti V ely for various values embedded flow parameter. From Table 

2, it is understood that, as increasing values of magnetic field parameter (M) results in considerable opposition to the flow 
in the form of a Lorenz drag which enhances the values of skin-friction coefficient, but there is a decrease in the wall 
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temperature gradient and the wall concentration gradient. The wall temperature gradient reduces as increase the values of 
radiation parameter R or dissipation Ec or heat generation parameter Q, while it is increases for increasing value of Prandtl 
number Pr. It is also observed that the increase in Schmidt number Sc or chemical reaction parameter Kr parameter lead to 
the increase in the dimensionless wall concentration gradient. 




Figure 2: Velocity Profiles for Varying Figure 3: Temperature Profiles for Varying 

Values of Magnetic Parameter (M) Values of Magnetic Parameter (M) 




Figure 4: Concentration Profiles for Varying Figure 5: Temperature Profiles for Varying 
Values of Magnetic Parameter (M) Values of Prandtl Parameter (Pr) 




Figure 6: Temperature Profiles for Varying Figure 7: Temperature Profiles for Varying 

Values of Radiation Parameter (R) Values of Viscous Dissipation Parameter (Ec) 
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Figure 8: Temperature Profiles for Varying Figure 9: Concentration Profiles for Varying) 
Values of Heat Generation Parameter (Q) Values of Schmidt Parameter (Sc) 




Figure 10: Concentration Profiles for Varying Reaction Rate Parameter (Kr) 



Table 1: A Comparison of Skin-Friction Coefficient •> v ' for Different 
Values of M for Fixed Values of Pr = R = Q = Ec = Sc = Kr = 0 





-/'(0) 


M 


Kameswaran et at. [19] 


Present 


0.0 


1.28181 


1.29038 


1.0 


1.62918 


1.63038 


2.0 


1.91262 


1.91285 


3.0 


2.15874 


2.15879 


4.0 


2.37937 


2.37938 



Table 2: Computation Showing f ^ , ® ^ and ^ ^ for 
Different Embedded Flow Parameter Values 



M 


Pr 


R 


Ec 


Q 


Sc 


Kr 


-/'(0) 


-0-(O) 


-^'(0) 


1.0 


0.71 


0.5 


0.1 


0.1 


0.6 


0.5 


1.63038 


0.518152 


0.898172 


2.0 


0.71 


0.5 


0.1 


0.1 


0.6 


0.5 


1.91285 


0.47744 


0.873141 


3.0 


0.71 


0.5 


0.1 


0.1 


0.6 


0.5 


2.15879 


0.445763 


0.854063 


1.0 


1.0 


0.5 


0.1 


0.1 


0.6 


0.5 


1.63038 


0.591577 


0.898172 


1.0 


2.0 


0.5 


0.1 


0.1 


0.6 


0.5 


1.63038 


0.832569 


0.898172 


1.0 


0.71 


1.0 


0.1 


0.1 


0.6 


0.5 


1.63038 


0.465892 


0.898172 
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1.0 


0.71 


2.0 


0.1 


0.1 


0.6 


0.5 


1.63038 


0.417967 


0.898172 


1.0 


0.71 


0.5 


0.5 


0.1 


0.6 


0.5 


1.63038 


0.355639 


0.898172 


1.0 


0.71 


0.5 


1.0 


0.1 


0.6 


0.5 


1.63038 


0.152498 


0.898172 


1.0 


0.71 


0.5 


0.1 


0.5 


0.6 


0.5 


1.63038 


0.354218 


0.898172 


1.0 


0.71 


0.5 


0.1 


1.0 


0.6 


0.5 


1.63038 


0.0846552 


0.898172 


1.0 


0.71 


0.5 


0.1 


0.1 


0.78 


0.5 


1.63038 


0.518152 


1.03456 


1.0 


0.71 


0.5 


0.1 


0.1 


1.0 


0.5 


1.63038 


0.518152 


1.18749 


1.0 


0.71 


0.5 


0.1 


0.1 


0.6 


1.0 


1.63038 


0.518152 


1.06315 


1.0 


0.71 


0.5 


0.1 


0.1 


0.6 


2.0 


1.63038 


0.518152 


1.33022 



CONCLUSIONS 

In this article we have studied the effects of radiation and viscous dissipation on heat and mass transfer from an 
exponentially stretching surface in the presence of heat generation and chemical reaction. The governing equations were 
solved numerically using the Runge-Kutta fourth order along shooting method. This has been shown to give accurate 
results. The effects of various physical parameters on the fluid properties, the skin-friction coefficient and the heat and 
mass transfer rates have been determined. We found that the effect of the magnetic parameter is to reduce the velocity of 
the fluid in the boundary layer region. It was also observed that the increase in values of M results in thickening of the 
species boundary layer. The chemical concentration boundary layer was found to decrease near the boundary with 
increasing reaction rate parameter and the Schmidt parameters. The heat transfer rates increases with an increasing of 
individual effects of the magnetic parameter M, the radiation parameter R, heat generation parameter Q and Eckert number 
Ec. The mass transfer rate increases with an increasing of Schmidt parameter Sc or the rate reaction parameter Kr. Also the 
numerical results obtained are agrees with previously reported case available in the literature Kameswaran et al. [27]. 
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